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ABSTRACT 

The projected virial theorem is applied to the case of multiple stellar populations in the nearby dwarf 
spheroidal galaxies. As each population must reside in the same gravitational potential, this provides strong 
constraints on the nature of the dark matter halo. We derive necessary conditions for two populations with 
Plummer or exponential surface brightnesses to reside in a cusped Navarro-Frenk- White (NFW) halo. We ap- 
ply our methods to the Sculptor dwarf spheroidal, and show that there is no NFW halo compatible with the 
energetics of the two populations. The dark halo must possess a core radius of ~ 120 pc for the virial solu- 
tions for the two populations to be consistent. This conclusion remains true, even if the effects of flattening or 
self-gravity of the stellar populations are included. 

Subject headings: Galaxies: dwarf - Galaxies: individual: Sculptor- Galaxies: kinematics and dynamics - dark 
matter 
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1. INTRODUCTION 

Using high resolution numerical simulations of structure 
formation, Navarro, Frenk & White (1995) first proposed that 
dark haloes have the universal form 



P(r) 



Po 



(r/r s )(l+r/r s ) 2 



(1) 



where r s is the scale-radius and po is the density normalisa- 
tion. Sometimes it is useful to define a halo radius r^, which 
is the radius within which the mean density is a factor of 
Ah = 200 times greater than the critical density. Then, the 
enclosed mass of the NFW profile is 



M(r) = Anp Q ti 



log(l + cx) 



1 + cx 



(2) 



where x = r/r^ and c = r^/r^. The two parameters of the 
NFW halo can be alternatively chosen as the concentration 
c and the mass M — M(rh). The NFW model remains the 
best two-parameter fit to simulation data, though Navarro et 
al. (2004) provided evidence that the three-parameter Einasto 
profile gives a still better match. 

Despite an extensive body of research, it is unclear whether 
the shapes of dark haloes inferred from observations are con- 
sistent with cosmological predictions. In galaxies like the 
Milky Way and M31, the dark halo profiles have been signifi- 
cantly affected by the cooling of baryons, which makes them 
unattractive candidates for testing the prediction. Therefore, 
research has concentrated on two very different classes of ob- 
jects, namely galaxy clusters and dwarf galaxies. In clusters, 
the high temperature and low density of the intra-cluster gas 
prevents efficient cooling, and so the baryons do not signifi- 
cantly affect the potential of the dark matter. There is reason- 
ably good evidence from cluster lensing that the profiles fol- 
low the NFW law, although with higher concentrations than 
predicted (e.g., Broadhurst et al. 2008). 

In dwarf galaxies at the present epoch, the baryonic com- 
ponent is feeble in comparison with the dark matter. Low sur- 
face brightness galaxies with HI and Ha rotation curves have 
been studied extensively by de Blok & Bosma (2002), who 
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conclude that cored rather than cusped haloes are favoured. 
The dwarf spheroidal galaxies have no gas, and are pressure- 
supported rather than rotation-supported, so are more difficult 
to assess. It was quickly realised that the radial velocities of 
thousands of bright giant stars could be gathered with mod- 
ern multi-object spectrographs (Kleyna et al. 2002; Walker et 
al. 2007). Dynamical modelling with either the Jeans equa- 
tions or phase space distribution functions then yields con- 
straints on the dark matter profile. However, this line of in- 
quiry seemed to have stalled with the realization of the degen- 
eracies inherent in the Jeans equations, making the datasets 
compatible with both cored and cusped dark halo models de- 
pending on the underlying assumptions (Evans et al. 2009; 
Walker et al. 2009). 

The problem has received new impetus with the growing 
realization that dwarf spheroidals often contain two or more 
stellar populations (Battaglia et al. 2006; Koch et al. 2006). 
This significantly strengthens any constraints from dynami- 
cal modelling, as now each sub-population must separately be 
in equilibrium in the same potential. There have been three 
recent and thought-provoking results. First, Battaglia et al. 
(2008) used the Jeans equations to model the velocity disper- 
sion profiles of the two populations in the Sculptor dSph, find- 
ing that a cored halo is preferred to an NFW halo, though the 
latter is still statistically consistent with the observations. Sec- 
ond, using the same dataset but more rigorous methods based 
on distribution functions, Amorisco & Evans (2012a) showed 
that NFW halos are rejected at high significance, performing 
substantially poorer than cored models. Finally, Walker & 
Penarrubia (2011) used the half-light radii of the two popu- 
lations in the Fornax and Sculptor dSphs to estimate the en- 
closed mass. This exploited earlier results showing that the 
mass within the half-light radius is robust against changes in 
the anisotropy (Walker et al. 2009; Wolf et al. 2010). Given 
the enclosed mass at two radii, the slope of density is inferred 
and shown to be in contradiction with NFW haloes. 

Here, we provide a simple line of reasoning based on the 
energies of the subpopulations. This brings the earlier argu- 
ments of Battaglia et al. (2008), Walker & Penarrubia (201 1) 
and Amorisco & Evans (2012a) into sharp focus. We show 
that the energetics of the subpopulations in one of the best 
studied dSphs, Sculptor, cannot be made consistent with an 
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NFW halo. 

2. THE PROJECTED VIRIAL THEOREM 

Let z be the line of sight direction. Then, the projected virial 
theorem relates the projected component of the pressure and 
potential energy tensors by the equation 2K\ 0S + W\ os = 0, 
where 



- t tS 



w l0 



-I 



vzd z (£>d 3 x, 



(3) 

Here, v is the luminosity density of a population with stellar 
mass-to-light ration Y* moving in the gravity field <Jx Un- 
der the assumption of spherical symmetry, the tensors be- 
come (e.g., Merrifield & Kent 1990) 
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(4) 
(5) 



where M(r) is the mass enclosed within radius r. These for- 
mulae involve the line of sight velocity second moment (v 2 os ) 
and the surface brightness //, which are directly accessible to 
observation. Notice that the velocity anisotropy of the stellar 
population does not enter, which gives the virial equations an 
immediate advantage over the Jeans equations. Notice also 
that virial quantities, which are gross volume integrals, can be 
computed more robustly than the gradients of observables re- 
quired by the Jeans equations. Furthermore, there is no prob- 
lem in dealing with rotation, as the integral /q os requires just 
the sum of square velocities in ordered and random motions 
along the line of sight. 

The only unknown in the projected virial equation is the to- 
tal gravitational potential <5(r) in which the population moves, 
or equivalently the enclosed mass M{r). If preferred, Wi os can 
be recast in terms of the total mass density p(r) which gener- 
ates d>, namely 
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Here, we have eliminated the luminosity density in terms of 
the surface brightness, which is an observable. 

The projected virial equation by itself is of course insuffi- 
cient to determine the potential or mass distribution directly. 
Suppose now there are two populations in equilibrium in the 
same total potential (e.g., Agnello & Evans 2012). In the 
dSphs, a younger metal-rich and an older metal-poor popu- 
lation with associated velocity dispersion profiles <x r and cr p 
and surface brightnesses // r and /i p (with corresponding lumi- 
nosity densities v r and v p ) are often present. For example, 
there is good evidence that Sculptor (Battaglia et al. 2008; 
Walker & Penarrubia 2011), Carina (Bono et al. 2010) and 
Fornax (Amorisco & Evans 2012b) contain at least two dis- 
tinct populations. Then, the projected virial theorem must be 
satisfied for each population, so we have 
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The only unknown is the potential, or equivalently, the en- 
closed mass M(r). This result is powerful enough to rule out 
some dark matter haloes. The physical reason is that a stellar 



population with velocity dispersion (v 2 os > offers the best con- 
straint on the potential at radii given by <D(r) « (v 2 os ). With 
more than one populations, the potential is restricted at more 
than one location. 

3. APPLICATION TO SCULPTOR 

3.1. Simple Results 

If the potential is heavily dominated by dark matter, we can 
neglect the self-gravity of the luminous components in a first 
approximation (relaxed in the next sub-section). In this case, 
both the kinetic energy tensor K\ os and the potential energy 
tensor Wi os scale linearly with the stellar mass-to-light ratio, 
which can then be eliminated. 

The surface brightness profiles of dSphs are well-fit by 
Plummer laws 



(1 + R 2 /R 2 h ) 



2\2 



(8) 



This assumption is not critical to our argument, but it does 
provide a useful illustrative model. We shall in any case 
sketch the extension to another commonly used profile, 
namely the exponential law. Fitting the Plummer law to 
the profiles of the metal-rich and metal-poor populations in 
Sculptor gives R^ = 230 ± 10 pc for the metal-rich and and 
Rh, P = 350 ± 10 pc for the metal-poor (consistent with the val- 
ues in Battaglia et al. (2008) and Amorisco & Evans (2012a)). 
A convenient fitting formula for the velocity dispersion profile 
is 



<vL> 



(1 +R 2 /a 2 ) ' 



(9) 



which gives a flattish profile with amplitude <tq out to a char- 
acteristic lengthscale a. This is typical of the velocity dis- 
persion profiles in, for example, Walker et al. (2007) or 
Battaglia et al. (2008). Again, this assumption is not nec- 
essary for our argument, as the relevant virial integrals can 
always be computed numerically. For Sculptor, fits provide 
(To , r = 8.7 + 1 .0 kms" 1 and a v = 240 ± 55 pc for the metal-rich 
and 0- o ,p = 10.9 ± 0.8 kms" 1 and a p = 1920 ± 850 pc for the 
metal-poor. 

With these laws, the kinetic energy integral is analytic with 
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The metal-rich limit is given by a r 
is given by a p » Rh, p . This gives 



^h,r, the metal-poor limit 
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The projected potential energy of Plummer profiles embedded 
in NFW halos is also exact, as 



W los = - 
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with A 2 = R\ + r 2 . With the values of Rh listed above, the 
virial ratio Wi os>p /Wi OS;r as a function of r s is increasing and 
bound from above. In particular, we have that 
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Figure 1. Left: Virial stripes for the two stellar populations in Sculptor in a cusped NFW potential, including the self-gravity of the stellar 
populations (Y* = 8). Purple shows the metal-rich population, blue the metal-poor population. In each stripe, the central line is the mean value 
of log 10 (p ), whilst the median and outer lines follow the l<x and 2cr deviations. Center and Right: Virial stripes for the two stellar populations 
in a NFW potential with small core, without (Y* = 0) and with (Y* = 8) self-gravity. 
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Table 1 

Minimum r s for two-sigma overlapping of the virial stripes. 



Hence, a necessary condition for a NFW halo to support two 
stellar populations with Plummer profiles is 



> 2 



^h,r 



(14) 



This is identical to eq (22) of Amorisco & Evans (2012a), 
derived under different assumptions. If, instead of Plummer 
profiles, exponential laws are used to fit the surface brightness 
profiles, then the numerical factor becomes 1 .9 instead of 2 in 
eqn (11). The analogue of eqn (13) is unchanged, so that the 
necessary condition for an NFW halo to support two stellar 
populations with exponential surface brightness profiles is 



CQ.i- 



> 1.9 



h. P 



(15) 



Using the best-fitting values provided above for Sculptor, it is 
immediate to check that the NFW potential is ruled out. Note 
that the constraints are simply the requirement that there is a 
NFW model with r s < oo. This is a much looser constraint 
than requiring consistency with an NFW model with a con- 
centration c m 20, as predicted by cold dark matter theories. 

3.2. The Virial Stripes 

The simple results already suggest that the energetics of the 
two populations are inconsistent with an NFW profile. How- 
ever, it is prudent to confirm this result numerically, discard- 
ing some of the simplifying assumptions made above. 

Since the measured profiles come with errors, we operate 
in the following manner. For each value of r s , we compute po 
separately for the two populations for many different photo- 
metric (8) and kinematic (9) profiles. We weight each result 
with the likelihood of the fit. Then, varying r s produces a 
virial stripe for each population in the (po, r s )-plane. If the 
two stripes overlap at 2<x at a particular r s , then the model for 
the potential is plausible at the 2<x level. Nothing prevents us 
from including the contribution of the luminous tracers to the 
potential as well. The virial equations then depend also on 



the stellar mass-to-light ratio T*, which may be different for 
the two populations. The projected potential energy Wi os has 
a contribution Wd m from the dark component and a correction 
W se if from the two luminous ones. For Plummer profiles, we 
have for the z— th population: 

W se i f j = n 2 Gp ,iRhj J] T+jpojRl jW (R hJ /R hJ ) , (16) 



with 



w(x) = 



t 3 [(5x 2 + 3)K(1 - x 2 ) - (x 2 + 7)E(1 - x 2 )] 
3(1 -x 2 ) 3 



(17) 



where K, E are complete elliptic integrals and Y*j is the lu- 
minous mass-to-light ratio of the j— th population. As the p ( y 
are given by number counts and not directly by luminosities, 
a common rescaling is applied to both populations such that 
the total luminosity is fixed at the observed value (taken from 
table 6 in Irwin & Hatzidimitriou 1995). 

The leftmost panel of Figure 1 shows the virial stripes for 
the two populations in Sculptor, excluding and including the 
effects of the self-gravity for the luminous component. The 
two stripes never overlap at the 2<x level. This confirms the 
result deduced from our simple argument in the previous sec- 
tion: there is no NFW halo compatible with the kinetic ener- 
gies of the two stellar components in Sculptor. The center and 
rightmost panels of Figure 1 show the virial stripes when the 
dark halo density is the simplest cored analogue of the NFW 
halo, namely 



cNFW = 



Po 



(6- 



2 /r 2 y' 2 (l+r 2 /r 2 ) 



(18) 



In this case, the stripes do overlap at the 2<x level provided the 
core radius r c = r s e is at least 150 pc, if the self-gravity of 
the stellar populations is neglected. Incorporating self-gravity 
causes the core radius to increase somewhat, as we see in 
Table 1. This shows how the minimum r s for 2<x overlap 
varies with changing Y* for different models. The first col- 
umn (Y* — > 0) stands for models in which the self-gravity of 
the stars is omitted. Since both half-light radii are smaller or 
equal to r s in the dark matter only case, adding self-gravity is 
expected to yield larger cores, as in fact is confirmed by the 
results in the Table. 

4. DISCUSSION AND CONCLUSIONS 

The arguments in this Letter show that the kinematics of 
multiple populations in dSphs provide a substantial challenge 
to the predictions of cold dark matter cosmogonies. In the 
case of one of the best studied dSphs, Sculptor, there is no 
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Figure 2. Virial stripes for the two stellar populations in Sculptor in 
a flattened NFW potential (q/g = 0.7 with Y* = 0). Purple shows 
the metal-rich population, blue the metal-poor population. In each 
stripe, the central line is the mean value of log 10 (po), whilst the me- 
dian and outer lines follow the lcr and 2a deviations. 

NFW dark halo that is compatible with the available photo- 
metric and kinematics data. The problem is very basic - the 
energetics of the metal-rich and metal-poor populations do not 
permit them to reside in the same NFW halo. Are there are 
any loop-holes? 

4.1. Reliability of Separation of Sub-Populations ? 

This Letter has used the separation proposed by Battaglia 
et al. (2008), who used a hard cut in metallicity to define 
the metal-poor ([Fe/H] < -1.7) and metal-rich populations 
([Fe/H] > -1.5). This is open to the objection that the 
kinematics and metallicity are interlocked, and so the cut 
may be subject to unperceived biases. The sub-populations 
can also be separated using a maximum likelihood approach 
with metallicity and kinematics treated jointly. We have 
checked our calculations using such a method on the sam- 
ple of stars in Tolstoy et al. (2004) (also used by Battaglia et 
al. 2008). Walker & Penarrubia (2011) have also tacked this 
problem, using a Bayesian likelihood estimator based on po- 
sitional, kinematical and chemical data to separate the sub- 
populations. Their algorithm also makes assumptions that 
may or may not be warranted (for example, the velocity dis- 
tributions are assumed isothermal). Nonetheless, although the 
details of the separation are different, Walker & Penarrubia 
(201 1) come to the same conclusion that the multiple popula- 
tions of Sculptor are inconsistent with an NFW halo. Different 
methods of separation, using different algorithms, come to the 
same answer, which offers comfort to those concerned about 
reliability. 

4.2. Steady state, well mixed populations ? 

An assumption in the virial theorem is that the stellar pop- 
ulations are in a steady-state and well-mixed. Recently, de 
Boer et al (2012) have provided a star-formation history for 
Sculptor. There was a peak in star formation ~ 13-14 Gyr 
ago, providing the bulk of the old, metal-poor population. Star 
formation continued till to ~ 7 Gyr ago. With the typical 
dark matter densities found from the virial stripes, the dynam- 
ical time-scale {AnGpa)~ l l 2 is at least one order of magnitude 
smaller than to. It seems that the two populations have in- 
deed reached a steady-state and that each of them is internally 
well-mixed. From the Hubble Space Telescope proper mo- 
tions, the orbital period of Sculptor is estimated as ~ 2.2 Gyr 
(Piatek et al. 2006). Its pericentric distance is ~ 68 kpc, and 



the time of the last pericentric passage was ~ 1 Gyr ago. As 
the dSph always remains in the outer parts of the Milky Way 
halo, it seems very improbable that tidal stirring can overturn 
the steady-state hypothesis. 

4.3. Flattening ? 

A drawback of all previous models of dSphs is that they are 
restricted to spherical symmetry. This is not the case for our 
virial arguments. 

The isophotes of the Sculptor dSph have a measured ellip- 
ticity of 0.3 (Irwin & Hatzidimitriou 1995). In fact, the flat- 
tening of the dark matter potential must be rounder than the 
flattening of the luminosity density, so spherical models are 
probably a good approximation. As a check, we generalize 
our virial arguments to flattened stellar populations (with axis 
ratio q) in a flattened NFW dark matter halo (with axis ratio g). 
Let us assume that Sculptor is viewed edge-on,and the density 
and kinematics of both stellar populations are stratified on sur- 
faces x 2 + y 2 q~ 2 - We also make the natural assumption that 
the minor axes of the dark matter and the luminous matter are 
aligned. Then, each kinetic energy term in eqn (7) scales by 
q, and so the ratio is left unchanged. Similarly, each potential 
energy term in eqn (7) is modulated by the same function of 
q/g in the integrand (c.f., the Flattening Theorem of Agnello 
& Evans 2012) and so the ratio is left largely unchanged. This 
suggests that our virial arguments will be robust against flat- 
tening. 

To confirm this, we carried out explicit calculations of the 
virial stripes for flattened models of Sculptor with q/g = 0.7, 
an example of which is shown in Fig 2. The virial stripes for 
the flattened NFW halo never overlap. Introducing a small 
core does permit the kinematics of the two populations to be 
consistent at the 2cr level providing the core is > 120 pc, a 
touch smaller than in the spherical case. 

We believe that the virial arguments are particularly robust 
and emphasise the important point that the energetics of the 
two populations in Sculptor are completely incompatible with 
al/r dark matter density cusp. For at least this dwarf galaxy, 
we can securely state that it either never had al/r density 
cusp, or if it did have one, feedback processes have now re- 
moved it. 

Of course, there has been much recent interest in converting 
cusped haloes into cored. A number of promising possibilities 
have been suggested, including impulsive mass loss from su- 
pernovae (Read & Gilmore 2005), winds and gas flows driven 
by supernovae (Mashchenko et al. 2006) and infalling bary- 
onic clumps (Cole et al. 2011). Which of these processes, if 
any, were responsible for the core in Sculptor is an important 
question for future work. 
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